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Abstract 

We derive and justify a normal form reduction of the nonlinear Schrodinger equation for a 
general pitchfork bifurcation of the symmetric bound state that occurs in a double-well symmetric 
potential. We prove persistence of normal form dynamics for both supercritical and subcritical 
pitchfork bifurcations in the time-dependent solutions of the nonlinear Schrodinger equation over 
long but finite time intervals. 



1 Introduction 

We consider the nonlinear Schrodinger (NLS) equation with a focusing power nonlinearity and an 
external potential (also known as the Gross-Pitaevskii equation), 

i^t = -'^xx + V{x)^-m^P^, (1.1) 

where ^{x,t) : M x M — )• C is the wave function, p G N is the nonlinearity power, and V{x) : R — )• M is 
the external, symmetric, double-well potential satisfying the following conditions: 

(HI) V{x) e L~(M) and xV'{x) G L°°(R); 

(H2) lim|,,|^^y(x) =0; 

(H3) V{-x) = V{x) for all x G M; 

(H4) Lq = —d"^ + V{x) has the lowest eigenvalue —E^ < 0; 

(H5) V{x) has a non-degenerate local maximum at x = and two minima at x = ±xq for some xq > 0. 
The easiest way to think about the double- well potential V{x) is to represent it with 

Vix) = ^[Vois-x) + Vois + x)], s>0, (1.2) 
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where the single- well potential Vq{x) satisfies (H1)-(H4) and has a global minimum at x = and no 
other extremum points. For sufficiently large s > s*, where s^, is the inflection point of Vq, that is, 
Vq{s^) = 0, the sum of two single-well potentials (jl.2|) becomes a double-well potential we would like 
to consider. 

The symmetric double-well potentials are used in the atomic physics of Bose-Einstein condensation 
PP through a combination of parabolic and periodic (optical lattice) potentials. Similar potentials were 
also examined in the context of nonlinear optics, e.g. in optically induced photo-refractive crystals [5] 
and in a structured annular core of an optical fiber [9] . Physical relevance and simplicity of the model 
make the topic fascinating for a mathematical research. 

Bifurcations of stationary states and their stability in the NLS equation (jl.ip under the assumptions 
(H1)~(H5) on the potential V{x) were recently considered by Kirr et al. [7]. 

Let = e^^^(j){x;E) be a stationary state such that (p{x; E) is a solution of the stationary 

nonlinear Schrodinger equation 

i-dl + V)^ - + E^ = 0. (1.3) 

Via standard regularity theory, if y S L°°(R), then any weak solution (f){-;E) £ H^{W) of the stationary 
equation (jl.3p belongs to ff^(R). Moreover, if -E ^ cr(Lo), then the solution (l){-;E) G H'^{R) decays 
exponentially fast to zero as |x| — t- oo. 

Existence of symmetric stationary states (p for any E > Eq bifurcating from the lowest eigenvalue 
— £^0 of the operator Lq = —d"^ + V{x) was first considered by Jeanjean and Stuart [3]. Kirr et al. [7] 
continued this research theme and obtained the following bifurcation theorem. 

Theorem 1 {Kirr et all. J^). Consider the stationary NLS equation lll.3\) with p > ^ and V{x) 
satisfying (H1)~(H5). 

(i) There exists a curve {Eq,oo) B E ^ <p{-;E) G //^(M) of positive symmetric states bifurcating 
from the zero solution at E = Eq. This curve undertakes the symmetry-breaking (pitchfork) bifurcation 
at a finite E^ G {Eq, oo), for which the second eigenvalue of the operator 

L^{E) = -dl + V{x) - {2p + l)<l?'P{x; E)+E (1.4) 

passes from positive values for E < E^ to negative values for E > E^. 

(a) Let 4>*{x) = (l){x;E^) be the positive symmetric state at the bifurcation point and G ff^(M) 
be the anti- symmetric eigenvector of L^{E^,) corresponding to the second eigenvalue X{E) such that 
X{E^) = and X'{E^) < 0. The curve {Eq,oo) 3 E ^ 4>{-'-,E) G ff^(M) intersects transversely at 
E = E^, with the curve of positive asymmetric states E i— t- ^±{-] E) G H^{M.) that extends to E > E^, 
if Q < and to E < E^ if Q> 0, where 

Q = 2p\2p + lf{^''r''^lLl\E,)ct>f-W.)L^ + \p{2p + l){2p-mlct?r^i^l)L2. (1.5) 
The asymmetric states (fj^ and tp^ are centered at the left and the right well of V , respectively. 

Orbital stability of the stationary state (^{x] E) in the NLS equation (II. ID depends on the number 
of negative eigenvalues of L^{E) and L^{E), where 

L.{E) = -dl + V{x)-(j)'^P{x;E) + E. (1.6) 
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Since L_{E)(p{E) = and (j){x; E) > for all x G M and E > Eq, the spectrum of L_[E) is non-negative 
for any E > Eq. This fact simplifies the stability analysis of the stationary states. 

Let us denote NsiE) = E)\\l^ and NaiE) = \\ip+{-, E)\\l^ = \\ip. {■; E)\\'j^2- In what follows, we 
always assume that 

N'g{E^) = 2{dE<f>*,(t>*)L^ > 0, where dE(t>*{x) = dE4>{x; E^), 

that is, Ns{E) is increasing near the bifurcation point E = E^. For example, this assumption is satisfied 
for large separation distance s between the two wells of V given by (jl.2p . because E^, ^ Eq as s ^ co 
and Ns{Eq) = 0. The following stability theorem was also proven in [7j. 

Theorem 2 [Kirr et all. Assume N'^{E^) > in addition to conditions of Theorem^ Then the 

symmetric state (p is orbitally stable for E < E^ and unstable for E > E^. If in addition, Q < 0, then 
Na{E) is an increasing function of E > E^, if S > and it is a decreasing function of E > E^ if S < 0, 
where 

S = K{E,) + (A'(i^.)||^*||i.)' • (1-7) 

Consequently, the asymmetric states (p± near E = E^ are orbitally stable for 5 > and unstable for 
5 < 0. 



For any potential V{x) represented by ()1.2p with sufficiently large s — )• oo, it was found in [7] that 
Q < for any p > \, hence, the stable symmetric state cj) for E < E^ becomes unstable for E > E^ 
and the asymmetric states (p± exist for E > E^. In the limit s — )■ oo, the boundary 5 = is equivalent 
to p = p^, where 

= ^ 3.3028. (1.8) 

li p < p^, the asymmetric states ip± are stable for E > E^. If p > p^, both symmetric and asymmetric 
states are unstable for E > E^. Therefore, we can classify the symmetry-breaking bifurcation at 
E = E^: as the supercritical (if 5 > 0) or the subcritical (if 5 < 0) pitchfork bifurcations with respect 
to the squared L^-norm, which is a conserved quantity of the NLS equation (|l.ip in time. The functions 
Ns{E) and Na{E) in the two different cases are shown schematically on Figure [H where stable branches 
are depicted by solid line and the unstable branches are depicted by dotted lines. 

The classification into the supercritical and subcritical pitchfork bifurcations is usually based on 
the analysis of the normal form equations obtained from the center manifold reductions and the near 
identity transformations. This analysis is the goal of this paper. We shall look at the long but finite 
temporal dynamics of the normal form equations, avoiding the complexity of the time evolution at 
infinite time intervals. To enable near identity transformations up to any polynomial order, we shall 
only consider the integer values of p. 

The normal form equations have been considered previously in a similar context. In the limit of large 
separation of the two potential wells, Kirr et al. derived a two-mode reduction of the NLS equation. 
Persistence of this reduction for periodic small-amplitude oscillations near stable stationary states was 
addressed by Marzuola & Weinstein [8] . Similar but more formal reduction to the two-mode equations 
was developed by Sacchetti [I^ using the semi-classical analysis. In comparison with [SlI^, Sacchetti 
|10j considered the defocusing version of the NLS equation, where the anti-symmetric stationary state 
undertakes a similar symmetry-breaking bifurcation. Based on the two-mode reduction, Sacchetti [TT] 
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Figure 1: Schematic representation of the supercritical (top) and subcritical (bottom) pitchfork bifur- 
cations in terms of the functions Ns{E) and Na{E). Unstable stationary states are shown by dotted 
curves. 
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also reported the same threshold as in (|1.8|) that separates the supercritical and subcritical pitchfork 
bifurcations. 

Unlike these previous works, we shall deal with a general symmetry-breaking bifurcation of the 
symmetric states. We develop a simple but robust analysis, which justifies a general normal form 
equation for the pitchfork bifurcation. Arbitrary bounded solutions of the normal form equation are 
proved to shadow dynamics of time-dependent solutions of the NLS equation (jl.ip near the stationary 
bound states for long but finite time intervals. Previously, only small-amplitude periodic solutions of the 
normal form equation were considered within the two-mode approximations in the large separation limit 
[8]. Also, compared with the sophisticated analysis based on Strichartz estimates and wave operators 
for the linear Schrodinger equations in [8], our analysis is only based on the spectral decompositions 
and Gronwall inequalities. Thus, we show how basic analytical methods can be used to treat time- 
dependent normal form equations for bifurcations in the nonlinear Schrodinger equations. Our main 
result is formulated in the following theorem, where we use notations of Theorems [J and [2j 

Theorem 3. Assume that there exists G {Eq^co) such that A'(£^*) < 0, N!,{E^) > 0, and Q < 0. 
Fix A/q and define AN = A/q — Ns{E^). There exists e > such that for any \A.N\ < e, there exists 
T > 0, G i/^ withAfo = \\'^o\\l2, and functions {0,E,A,B) G C^{[0,T];R'^) such that the NLS 
equation lll.l]) admits a solution ^ G C{[0,T]; (W)) with '${x,0) = "^o{^) ii^ the form 

^{x, t) = e^^W [(Pix; E{t)) + E{t)) + iBxix; E{t))] + t), 

where if) {x]E) andx{x]E) satisfy 

L+{E)i; = -A\E)x, L^{E)x = i^, 

subject to the normalization (X;^)l2 = 1? with the asymptotic expansion, 

A^{E) = ->!{E^)U4l2{E-E^) + 0{E-E^f as E ^ E^. 

Moreover, there are positive constants Cq, Ci, C2, C3, and C4 such that T < Co|AA^|~^/^, 

< C'llAiVl, \e{t)-E^\<C2\AN\, \E{t) - E^\ <C3\AN\, for ah t G [0, T], 

and the trajectories of {A, B) in the ellipsoidal domain, 

D = {{A,B) G : A^ + \AN\-^B'^ < Ci\AN\] , 

are homeomorphic to those of the second-order system, 

1^ = ^' (19) 
1 B = {-\'{E,)\\i;4l,{AN)A+ QSA^) /N',{E,). ^ ' > 

Trajectories of the second-order system (jl.9p with A'(i?*) < and N'^^E^) > are shown on Figure 
[2] for four distinct cases of different values of AA^ and QS. 

The article is organized as follows. In Section 2, we use Theorem[Ui) and derive modulation equations 
for dynamics of time-dependent solutions of the NLS equation near the stationary bound states at 
the onset of the symmetry-breaking (pitchfork) bifurcation. In Section 3, we consider the stationary 
modulation equations and recover the results of Theorems HJii) and [2] from our system of equations. 
In Section 4, we justify the dynamics of the time-dependent modulation equations and give a proof of 
Theorem [3l 
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Figure 2: Trajectories of the second-order system (II. 9p on the phase plane {A, B) for AA^ > and 
Q5 < (top left); A7V > and Q5 > (top right); AiV < and Q5 < (bottom left); AiV < and 
Q5 < (bottom right). 
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2 Modulation equations for dynamics of bound states 



We shall derive a set of modulation equations which describe temporal dynamics of solutions of the 
NLS equation (jl.ip near the stationary bound states at the onset of the symmetry-breaking bifurcation. 
We only use the statement of Theorem [T][i) on the existence of the symmetry-breaking bifurcation for 
the symmetric stationary state of the NLS equation (jl.ip under assumptions (H1)-(H5) and p G N. 

2.1 Primary decomposition near the symmetric stationary state 

Let (j){x; E) be a solution of the stationary NLS equation ()1.3p with properties 

(^(•; E) G i?2(M) : (^(-x; E) = (p{x; E) > for aU x G M. 

It is stated in Theoremd^i) that the curve E (pi-; E) G H'^{W) exist for all E G {Eq, oo). If p G N, 
this curve is actually C°° by the bootstrapping arguments. 

We shall consider a solution of the NLS equation (jl.lj) in the form 

^(x, t) = e*^W [0(x; + u{x, t) + iw{x, t)] , (2.1) 

where {E, 9) are coordinates of the stationary state. Direct substitution of (|2.ip into (jl.ip shows that 
the real functions (n, w) satisfy the system of time evolution equations 

ut = L^w + N^{u,w) + {e-E)w-EdE(l), (2.2) 
-wt = L+u + N+{u,w) + {e- E){(l) + u), (2.3) 

where L_|_ and L_ are defined by ()1.4p and ()1.6p and the nonlinear terms are given explicitly by 

N_{u,w) = -w;[(02 + 2<^n + 'u2 + u;2)P-</,2p]. 
For any p G N, we can use the Taylor series expansions 
N+{u,w) = -p{2p+l)(t)'^'P"^u^ -p(j?P~^w^ 

-\p{2p + l){2p - l)(/)2P-2^3 _ p^2p - l)02p-2^^2 ^ 0(^2 ^ ^2-^2^ ^^A) 

O 

N_{u,w) = -2p(t)^P~^uw-p{2p-l)(l?'P"'^u^w-p4)^^~'^w^ + 0(u^ + w'^f. (2.5) 

To determine {E,0) uniquely in the neighborhood of the stationary state (for small u and w), we 
add the standard conditions of symplectic orthogonality 

{cl),u) L2=0, {dE(p,w)L2=0, (2.6) 

where we recall that 

L_(l> = 0, L+dE(l) = -((>■ (2.7) 
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Under symplectic orthogonality conditions (j2.6|) . the rate of changes of {E, 6) are uniquely determined 
from the projection equations 



{dE<f),<i)- u)l2 -{(j),w)L2 

-{dl4>,w)L2 {dE4>,4> + u)L2 

We shall now study eigenvectors at the onset of the symmetry-breaking bifurcation in order to build 
a frame for the secondary decomposition of the perturbations (n, w) near these eigenvectors. 



E 
9- E 



{4>,N_{u,w))l2 
-{dE(/),N+{u,w))L2 



2.2 Linear eigenvectors 

It is stated in Theorem [TJi) that there exists a bifurcation value E^ G {Eq,oo) such that the second 
eigenvalue X{E) of L^{E) satisfies A(£'*) = 0. We shall denote (l)*{x) = at the bifurcation 

value E = E^:. In many cases, we will suppress the x-argument in the function (j)(x;E) to underline 
the E'-dependence of this function. In this setting, we have the following result. 

Lemma 1. There exist odd functions tp*,x* ^ -ff^(M) such that 

L+{E,)^, = 0, L.{E,)x* = A, (x*,V'*)l2 = 1. (2.9) 

Moreover, 

X'{E.) = 1 - 2pi2p + (2.10) 

Proof. Let g(x;E) be an eigenfunction of L^{E) with eigenvalue X{E). By Sturm's Theorem, g{x;E) 
is odd in x because X{E) is the second eigenvalue of L+{E). Let ip^{x) = g{x;E^). So, is an odd 
function such that Lj^[E^,)ip^: = A(£^^,)V'* = 0. 

Since (t)*{x) = (j){x] £'*) is even, positive and L-{E^,)(j)^ = 0, we see that zero is the lowest eigenvalue 
of L-{Etf). Therefore, it is a simple eigenvalue. Because is odd, there is an odd function x* G -ff^(M) 
such that L^{E^)x* = tp*- On the other hand, we have 

(x*,^*)l2 = {LZ^{E^)tp*,'il)*)L^ > 0. 

By rescaling and x*, we get {x*,'>P*)l2 = 1. This completes the proof of the first part of the lemma. 

To prove (I2.10p . we note that due to the smooth continuation of (p{x; E) across E'*, we can compute 
explicitly, 

L\{E,) = 1 - 2p(2p + l)(fiP-^dE(l}., L'^{E,) = 1 - 2p(t)lP-^dE(p.. (2.11) 
By differentiating the relation L^{E)g{E) = \{E)g{E) at E = we get 

L'^{E^)ij^ + L+{E^)dEi'* = \'{E^)^^. 

Taking the inner product of this equation with ip^, we get 



L2 



y*||^2 II r* 112^2 

This completes the proof of the lemma. □ 
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We would like now to extend the functions {ip*,x*) ^ the eigenvectors of the linearized system 
associated with the time-dependent system (j2.2p and (j2.3p near E = E^. Note that the eigenvectors of 
the linearized (non-self-adjoint) system are different from the eigenvector g of the (self-adjoint) operator 
Lj^{E) introduced in the proof of Lemma [TJ The following lemma gives the extension of {^tp-t,^ x*) near 
E = Ej,. 

Lemma 2. There exists sufficiently small e > such that for all \E — E^\ < e, there exists a small 
eigenvalue A{E) of the linearized system 

L+{E)i^{E) = -K\E)x{E), L^{E)x{E) = ^{E), (2.12) 

where the eigenvector and eigenvalue satisfy the asymptotic expansion, 

i^{E)=ij, + OH2{E-E,), x{E) = x* + Oh2{E-E,), (x(^), V(^))l2 = 1, (2-13) 

and 

A^{E) = -X'{E,)\\i^4l2{E - E,) + 0{E - E,f . (2.14) 
Consequently, if \'{E^) < 0, the eigenvalue A{E) is real for E > E^ and purely imaginary for E < E^,. 

Proof. Recall that L-{E:^)(j)^ = 0, L^{E:^)7p^ = and L-{E^)x* = ip*, where both ip^ and x* odd 
in X and (p^ is even in x. By the perturbation theory for isolated eigenvalues, there exists a solution 
for A{E), 'il^{E), and x{E) in the system (|2.12p such that the eigenvectors satisfy the expansion (j2.13p . 

By taking the inner product of L+{E)ij{E) = -A^{E)x{E) with ip{E) and using (|2l0]l . (fZTT]) . and 
(j2.13p . we get (I2.14|) . So, it follows that if \'{E^) < the eigenvalue A(-E') is real for E > E^ and purely 
imaginary for E < E^,. 

On the other hand, for small \E — E^\, we have 

{x{E),^P{E))l2 = {x*,^P*)l^ + O{E-E,) = 1 + O{E-E,)>0. 

We can hence normalize ip{E) and xi^) such that {x{E),7p{E)) ^2 = 1. This completes the proof of 
the lemma. □ 

Remark 1. Under the normalization (x(-E), ^(-E))i2 = 1, the L? norms of ijj and x o,re no longer 
normalized to unity, in comparison with the normalization used in 

2.3 Secondary decomposition near the linear eigenvectors 

Let us now decompose the perturbation terms into 

u{x,t) = A{t)i){x;E) + U{x,t), w{x,t) = B{t)x{x]E) + W{x,t), (2.15) 

where {A, B) are coordinates of the decomposition and {U, W) are the remainder terms. The linear 
eigenvectors {ip,x) ^^e solutions of the linearized system (|2.12p for E near E^. The remainder terms 
([/, W) are required to satisfy the conditions of symplectic orthogonality 

(0, [7)^2=0, {dE<t),W)L2=^, {x,U)l2=0, {iP,W)l^=0. (2.16) 
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Substitution of (|2.15|) into (|2.2p - (|2.c{|) show that ([/, W) satisfy the time-evolution equations, 



Ut = L_W + N_{A4: + U,Bx + W) 

+ {e- E){Bx + W)- E{dE(t> + Adsi^) -{A- B)i^, 
-Wt = L^U + N+{A^ + U,Bx + W) 

+ {e-E){<p + Ai) + U) + BEOeX +{b- h^A)x. 



(2.17) 
(2.18) 



Under the orthogonahty conditions ()2.16p . the rate of changes of {E, 6, A, B) are uniquely determined 
from the projection equations 



M 



E 




{(t>,N^{A^ + U,Bx + W))L^ 


e-E 




-{dE(t^,N+{Ai, + U,Bx + W))L2 



where 



and 



M 



-{dlcP,W)L2 {dEcPA + U)L^ 



(2.19) 
(2.20) 



A-B 
B-K^A 



{x,N^{A^ + U,Bx + W))l2 

+E{{dEX, U)l2 - A{dE^,x)L^) + {e- E){B\\x\\l2 + {x,W)l2) 
-{i^,N+{A^ + U,Bx + W))L^ 

+E{{dE^I^,W)L2 - B{dEX,i^)L^) - {e-E){A\\iP\\l2 + (^,c/)l2; 



(2.21) 
(2.22) 



In Section 4, we shall control the dynamics of small {U,W), [E — Ejf,6 — E^,), and {A,B) in the 
system (j2.17p - (j2.22p on long but finite time intervals. 



2.4 Conserved quantities 

The NLS equation (jl.ip admits two conserved quantities given by 

^f[i^l] = I \qi{x,t)\'^dx 

Jr 



and 



\^,{x,t)\^ + V{x)\^{x,t)\' 



1 



\^{x,t)\ 



2p+2 



dx. 



(2.23) 



(2.24) 



P + l 

They are referred to as the power and the Hamiltonian respectively. 

Let Ns{E) = 7\A[(^(-; E)] and Hs{E) = [(/)(•; E)]. If ^'o is an initial datum for the solution ^' of the 
NLS equation (jl.ip . we define 

AAo = M[^q] and Ho = ?^[*o]- 
Substitution of and (f2T5]) into ^OJM and (12:21) gives 



AAo = Ns{E) + / [{Ai) + Uf + {Bx + Wf] dx 
Jr 
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and 



= Hs{E)+ f [{AiP^ + U^f + {Bxx + W^f + V{A^ + Uf + V{Bx + Wf]dx 



P + 1 Jr L 



j [((0 + Aij + Uf + {Bx + WfY^^ - (/>2p+2 _ 2{p + + U)] dx, 



where we have used the stationary equation (jl.3p for (j) and the symplectic orthogonahty (|2.16p . 
Remark 2. By direct computation, we can verify that 



for any E, for which (p{-;E) G H'^{]^) exists. 

3 Stationary normal-form equation 

We shall recover the results of Theorems dl^ii) and [2] on the existence and stability of stationary states 
from the system of time evolution equations (I2.17p ~ (l2.22p . Theorems [1] and [2] were originally proved 
in [7] with the Lyapunov-Schmidt decomposition method that relies on an orthogonal decomposition 
with respect to the self-adjoint operator L^{E). On the other hand, the decomposition used in the 
derivation of system (|2.17p - (j2.22p relies on the symplectic orthogonality conditions (j2.16p . Therefore, 
computations of this section provide an alternative proof of Theorems Hl^ii) and [2j 

3.1 Alternative proof of Theorems [T](ii) and [2] 

We start with the simplification of the system (j2.17p - ()2.22p for stationary solutions of the NLS equation 
(jl.ip . Because of the symplectic orthogonality conditions (I2.16p . we define the constrained H^-space, 



where the subscript indicates that the orthogonal projections are -E-dependent. The stationary solutions 
of the system (I2.17p - (l2.22p near E = E^ are described by the following theorem. 

Theorem 4. Assume that Ng(E^,) ^ 0. There exists sufficiently small e > such that for all \E — E:^\ < 
e, the nonlinear Schrddinger equation admits a stationary solution 



HiiE) + EN',iE) = 0, 



(2.25) 



Hl; = {Ue H'{R) : {<P{E), U)l2 = {x{E), U)l2 = O} 



(3.1) 



= e**^ [^{x; E) + Ai;{x- E) + U] 



(3.2) 



where 



while U € H]^ 




(3.3) 



L+{E)U = G{A, E, U), mE),G{A, E, U))l^ = 0. 



(3.4) 
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with 

GiA, E, U) = A^Ax + ^'"^/f't^ltt.f'^^^^^ l^ + Ai. + U)- N^iA^ + U,0). (3.5) 
Moreover, there exist positive constants Ci and C2 such that 

\\U\\h2 <Ci\E-E^\, A^<C2\E-E^\. (3.6) 



Proof. For real- valued stationary solutions, we can set B = and 14^ = in the system (|2.17|) - (|2.22p . 
which give w = and (Atl; + U,0) = 0. The modulation equations (j2.19p - (j2.22p become degenerate 
and can be rewritten in the form, 

{dE(t>A-U)L2E = 0, 
{dEct),cl) + U)L2{e-E) = -{dECI).N+{All, + U,Q))L2, 

A + E{-{dEX,U)L2+A{dEi^,x)L^) = 0, 

-K^A + {9-E){AU\\h + {^,U)L2) = -{iP,N+{Aij + UML2- 

Assuming Ni.{E^:) / and the smallness of ||C/||l2 for small \E — E^\, we get {dEfp^ 4' ~ U)l2 / for 
small \E — E^\. So, from the first equation, it follows that E = Q. Then, from the third equation, we 
infer that ^ = 0. The second equation gives (j3.3p with the correspondence 6 = 8. The fourth equation 
gives (V', G{A, E, U))i2 = 0, where G is defined by (13. Sp and the normalization x)l^ = 1 is used. 

The time-evolution system (I2.17p - ()2.18p implies that U becomes time-independent and satisfies the 
stationary equation L^U = G{A, E, U). So, the system ()3.4p is verified and we shall prove the existence 
and uniqueness of small solutions of this system satisfying the estimates (j3.6p for small | — | . 

Since Lj^{E^:)ip^: = and ~^ ^* L?'{^) as -E — )• i?*, the Implicit Function Theorem gives the 

existence of a unique local map 

3 {A,E)^U e H'^i^) near (^, ^) = (0, E,), (3.7) 

such that U satisfies equation L+{E)U = G{A,E,U) subject to the constraint {x{E),U)i2 = 0, 
provided that ('(/'(£'), G(^, ii^, J7)) ^2 =0. Let C/a.e denote this map for small A and | — | . The map 
is C°° if p G N. 

In addition, it follows from equation 

L+{E)dE(l>{E) = -(t>{E) (3.8) 

that Ua,e satisfies {4>{E),Ua,e) l'^ = under the constraint {dE(p{E),G{A, E,U)) ^2 = 0, which is 
identically satisfied. Therefore, Ua,e G H^, according to the definition (|3.ip . 

We note that G{A, E, U) is quadratic in A as j4 — )■ 0. Therefore, we proceed with a near identity 
transformation for the map (13. 7p . 

Ua,e = A^Q{x-E) + Oh2{A''), (3.9) 
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where G G is a unique solution of the inhomogeneous equation, 

L^{E)Q = p{2p + 1)02^- - Pi2p + ^) ^^''tfrlf^'^ ^- (3.10) 

It remains to control the value of A in terms of the small value of \E — E^:\. Substituting ()2.14p and 
dSJ]) to equation {iP{E),G{A,E,U))l2 = 0, we obtain 

X'{E^)\\i;4l2{E - E^)A - QA^ + 0{A\A^{E - E^), {E - E^f) = 0, (3.11) 

where 

Q = ^p{2p+l){2p-l){ijlAf-^4'l)L^+2p{2p + l){i^lcf'j>-^Q,)L2 

-p{2p + i) — , . h/'* L2- 

{dE(p*,<?*)L-2 

Therefore, either ^ = (and Ua,e = 0) or ^ is a non-zero root of equation (j3.1ip that satisfies the 
second estimate (|3.6p . Thanks to the expansion p.9p . the first estimate p.6p is also satisfied. This 
concludes the proof of the theorem. □ 

We will now show that the results following from Theorem H] and the normal form equation (|3.1ip 
are equivalent to the results of Theorems (U^ii) and [2j Using (|3.8p , we let 

G = ^+p(2p+ l)(^^M^^^a^<A, 9 = p{2p+l)L-\E)<t>''^-',l^\ 
On the other hand, expanding (|3.3p gives 

S = E + p{2p + i)^2 (gi.y^^- V^)L2 ^ 0(^3)_ (3 

Using (|2.10p . we conclude that the normal- form equation (|3.1ip is equivalent to equation 

\'{E,)\\il:4l^{E - E,)A - QA^ + 0(A^ ^^(^ _ _ ^^^2^ ^ ^^ -^^^ 

where Q is given by (jl.Sp and is a renormalized parameter of the stationary state (f). 

We can see from (|3.13p that besides zero solution A = that corresponds to the symmetric state 
(l){x] £) with £ = E, there are two nonzero solutions that correspond to the asymmetric states, 

(p±{x;£) = (l){x; E) ± Aip{x] E) + A^e{x; E) + Oh2{A^), (3.14) 

where £ is related to E by the expansion (I3.12p and yl is a positive root of the stationary normal-form 
equation ()3.13p provided that sign((<5 — E^)Q) = —1 if A'(£'*) < 0. If Q < 0, the asymmetric states 
exist for <S > If Q > 0, the asymmetric states exist for £ < E^. This concludes the alternative 
proof of Theorem [Hii) . 

To establish the alternative proof of Theorem [2l we need the following result. 
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Lemma 3. Let ip± he the asymmetric states JXT^P ^hat exist for sign {{£■ — E^)Q) = —1 near £ = 
and define 

L+{A) = -dl + V{x) - {2p + + £{A), 

where £{A) is given by the expansion ^3.13\) . Then, the second eigenvalue of L^{A) is positive for 
Q < and negative for Q > 0. 

Proof. Using (j3.1ip and (j3.14p . we expand ip± for small A by 

^± = ^,± A^, + A^ (q, + ^ dE(p^i + [A"). (3.15) 

Let h{A) be the eigenfunction of -^^+(^4) for the second eigenvalue ^^{A). We have /i(0) = ip^, and 
^(0) = is a simple eigenvalue, so that the analytic perturbation theory for self-adjoint operators 
applies. Using (I3.13P and (I3.15p . we compute 

L'+(0) = -2p(2p+l),^2p-i^^ 

and 



L'I(O) = -2p{2p + l){2p-l)<t>lP-^i,l-Ap{2p + l)ct>f-^i,, [q, + ^^^—^^j^S^,^,^ + „2 



^2 / A'(£^*;||</v*||^2 

Algorithmic computations yield /i'(0) = 26^, m'(0) = 0, and, after tedious computations, 

(L'|(0)V^*+4LV(0)^*,V'*)l2 _ 4Q 



1^2 IIV'*ll^2 



Therefore, ii{A) > for smah yl if Q < and n{A) < for smah A if Q > 0. □ 

Let Ns{E) = \\(j){-;E)\\\^ and Na{£) = \\^p±{-;£)\\\2- Using equations (fXT2]l and (|3l^ . we expand 
Na{£) in the power series, 

Na{£) = Ns{E) + A^\mh+0{A^) 

= Ns{E,) + N'^{E,){E - E,) + A^\\iP4l2 + 0((S - E,f, {E - E,)A^, A^) 
= Ns{E,)+S{£ - E,) + 0{£ - E,f, 

where 5 is given by (|1.7p . Assuming that Q < 0, the asymmetric states exist for £ > E^. If 5 > 0, 
then Na{£) increases with £, whereas if 5 < 0, then Na{£) decreases with £. 

Orbital stability and instability of asymmetric stationary states follows from the classical theorem of 
Grillakis, Shatah, & Strauss [2] because Lemma [3] shows that the operator Lj^{A) linearized at (p±{£) 
has one negative eigenvalue if Q < and £ > E^. On the other hand, the symmetric state </>(£') is 
unstable for > i?* by a theorem of Grillakis [3] because Theorem [T][i) shows that L^(E) linearized 
at (p{E) has two negative eigenvalues for E > E^,. This concludes the alternative proof of Theorem [2j 
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3.2 Limit of large separation of potential wells 

The case of large separation of potential wells, when s — t- cxo in the double- well potential (|1.2p . gives 
a good example of explicit computations of numerical coefficients Q and S. Knowing these numerical 
coefficients enables the explicit classification of the stationary states in Theorem [1] and [2] and verifies 
the conditions of Theorem [3l The following theorem gives the asymptotic result when s — )• oo. 

Theorem 5. Let V be given by (II. 2|) . There exists sufficiently large sq > such that for all s > sq, 
we have N',{E^) > 0, X'{E^) <0, Q<0, and 

S > 0, if p < p^, and S < if p < p*, 

where is the positive root of the equation 1 + 3p — p^ = 0. 



Proof. We recall from [7] that 

Eo and t/;^ ^ C^^^ in L°°(M) as s — )• oo, 
where C* > is a normalization constant. Using the exact identity 

we can hence simplify the expression (jl.Sp to the form, 

Q^-^p(p+l)(2p + l)a'||<A*||f2p+2 as s^oo. (3.16) 

Therefore, Q < for any p € N if s is sufficiently large. 

Because E^, — )• Eq, we can approximate (p^, and dE4'* using the small-amplitude expansion for the 
symmetric states of the stationary equation (|1.3p . 

II i ||2p+2 

cPix; E) = aMx) + 0{a'+'n, E = Eo + + 0{a'n, 

\m\\L2 

where 0o G -ff^(M) is the eigenfunction of the operator Lq = —d"^ + V{x) for the lowest eigenvalue —Eq 
and a G M is a small parameter of the expansion. As a result, we obtain 

{dE<t^*,4'f^^)L^^^\Mh, as s^oo. (3.17) 
Ip 



and 



K{E,) = 2{dE(t>*A*)L^ ^ - „ ,*2L+2 ^ as S^OO. (3.18) 

In addition, it follows from (j2.10p and ()3.17p that 

A'(^=,) -2p as s ^ oo. (3.19) 
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Therefore, N!,{E^) > and A'(£;*) < as s ^ oo. 

Substituting (|3.16p - (j3.19p into the expression (|1.7p . we obtain 



(l + 3p-/)||(^,r 
p(l+p)(l + 2p)||0,||f++2 



as s — )• DO. (3.20) 



Therefore, 5 > for p < and S < for p > p^, where is given by the positive root of l+3p—p^ = 0, 
that is, by ([L8])- □ 



4 Time-dependent normal form equations 

We rewrite in the abstract form the modulation equations (I2.19p - ()2.22p for {9, E, A, B), 

( 9 - E = Rg{E,A,B,U,W), 
E = REiE,A,B,U,W), 
A - B = Ra{E,A,B,U,W), 
^ B - A^{E)A = Rb{E,A,B,U,W), 

and the system (I2.17p - (j2.18p for the remainder terms {U,W), 

Ut = L^{E)W + Ru{E,A,B,U,W), 
-Wt = L+{E)U + Rw{E, A, B, U, W), 



(4.1) 



(4.2) 



where Rg, Re, Ra, Rb, Ru and Rw are some functionals on the solution. These functionals can be 
computed explicitly. Indeed, it follows from ()2.19p that 



Re 
Re 



{<t),N.{Ai; + U,Bx + W))L2 
-{dE(k,N+[Aij + U,Bx + W))l2 



(4.3) 



where matrix M is given by (I2.20|) . If ||[/||j;^2 , ||VF||2,2 ^ 1 and Ng{E^) / 0, then M invertible and 



= {dE<t^A)ll 



1 

1 



+ 0{\\U\\l2 + \\W 



Il2; 



(4.4) 



On the other hand, equations (j2.2ip and ()2.22p yield 

Ra = ix, N-{Ai; + U,Bx + W))^^ + ReHOeX, U)l^ - A{dE^P, x)l^) + 

Re{B\\x\\l2 + {x,W)L2), 

-{tlj,N+{Ai^ + U,Bx + W))L2 +RE{{dEi',W)L2 - B{dEX,i^)L^) 



Rb 



-ReiA\ml, + {^P,U)L2] 



Now, the system (j2.17p - (|2.18p simplified with the modulation equations (j4.ip yields 

i Ru = N^iA^P + U,Bx + W) + Re{Bx + W) - RE{dE(t> + ASe^^) - Ra^, 
I Rw = N+{A^ + U,Bx + W) + Rg{<p + Ai; + U) + BReOeX + RbX- 



(4.5) 
(4.6) 

(4.7) 
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Remark 3. If A = B = 0, then U = W = is an invariant solution of the system [J^, which give 
zero values of Re, Re, Ra, o,nd Rb in the system ^.1^ for any E. 



Observation of Remark [3] inspires us to consider the power series expansions for solutions of the 
systems (14. ip and (14. 2p . Taking into account the spatial symmetry of eigenfunctions, we can see that 
[/, R0, and Re are quadratic with respect to {A,B), whereas Ra and Rb are cubic with respect 
to {A,B). Therefore, we write 

' e-E = Ci{E)A'^ + C2{E)B'^ + Re{E, A, B, U, W), 

E = C3{E)AB + Re{E,A,B,U,W), 

A-B = C4iE)A'^B + C5iE)B'^ + RA{E,A,B,U,W), 
^ B - A^{E)A = Cg{E)A^ + C7{E)AB^ + Rb{E, A, B, U, W), 



\ U = A'^@{x; E) + 52A(x; E) + A^Ui{x- E) + AB'^lhix; E) + U{x, t), 

\ W = ABT{x;E) + A^BWi{x;E)+B^W2{x;E) + W{x,t), ^ ' ' 

where R0, Re, Ra, and Rb are the error terms, whereas U and W are the remainder terms. 

Remark 4. B = and W = yield the invariant reduction of the stationary system ^.8^ to the 
algebraic equation, 

-A\E)A = C6{E)A^ + Rb{E, A, 0, U, 0), 
from which it follows that Cq{E^,) = Q in notations of Section 3.1. 



Let us first explicitly compute the coefficients Ci , C2, . . . C7 and determine the functions G, A, . . . , VF2. 
We shall then estimate the error and remainder terms in (j4.8p and (j4.9p as quadric with respect to 
{A, B). Working in a small neighborhood of (0, 0) on the phase plane (A, B) and using | AA^| as a small 
parameter, we consider an ellipsoidal region on the (A, i3)-plane such that 

3C>0: A^ + \/^N\-^B^ <C\/\N\. (4.10) 



Let T > be the maximal time until which we consider solutions of the modulation equations (j4.8p 
in the domain (j4.10p . We assume (and prove in Section 4.4) that there are positive constants Cq, Ci, 
and C2 such that 

T < Co|AiV|~i/^ (4.11) 

and 

1^-^*1 < Ci|AiV|, 1^-^,1 < C72|A7V|. (4.12) 



The following theorem provides the control of the error terms of the system (j4.8p and the remainder 
terms of the decomposition 



Theorem 6. Assume ^AD^-^TW^. There exists e > such that for any |AA^| < e, there are positive 
constants Ci and C2 such that 



sup ( \\Ui;t)\\j,. + \\Wi;t)\\j,. ) < CiiANY (4.13) 

te[o,T] 
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and 

^2 



sup ( \Re\ + \Re\ + \Ra\ + \Rb\ ) < C2{ANy. (4.14) 

ig[0,T] 



The proof of Theorem [U] is given in Sections 4.1 and 4.2. 



4.1 Power series expansions 



For exphcit computations, we use the power series expansions (j2.4p - (j2.5p and the decompositions (|2.15p 
and (|4.9p to expand 



-p{2p + l),/.2f- V'^' - pc^^^-'x^B^ - p{2p + 1) (^^l-l^^P-^i;^ + 2(/.2P- 



- p {{2p - l)02p-Vx' + 2(2p + 1)</.2P- VA + 2</)2p-i;^r) AB^ + N+, 
N_ = -2p(l?P-^^xAB - p {^"^P-^x^ + 2(l)'^P-^xA) B^ 

- p {{2p - i)(i?p-'^xi^'^ + 2(/)^^" Vr + 2(l?P~^xA) A^B + iV_, 

where and are of the form 

iV+, iV„ = 0{{A^ + B^f + {A + B){U + W) + U^ + W^). (4.15) 
From equations (j4.7p . (j4.8p . and (j4.9p we have 

r = h,i{E)AB + f2,i{E)A^B + fo,3{E)B^ + Fu{E, A, B, U, W), 



Rw = 92,o(.E)A^ + go,2{E)B^ + g3,o(.E)A^ + 5i,2(i?)^i?' + Fw{E, A, B, U, W), 



where 



and 



52,0 = -p(2p+l)02p-V' + Ci(^, 

50,2 = -p</''^-'x' + 

/2,i = CiX - C^dE^P - C^^P - p{2p - 1)02P"2^^2 _ 2p^2p-i^r - 2p(t?P-\Q, 

51,2 = C2V' + C3ai,x + C7X-p((2p-l)</''P-W + 2(2p + l)02p-VA + 2,/.2p-V) 

53,0 = CiV' + - P(2p + 1) (\{2p - l)02p-2^3 ^ 202P" Ve 



Fc/ = (C7i^2^C25^)VF + (Sx + W^)i?e- (ai=;0 + Aai5^)i2i5-i?AV' 

-2p</.2p-i (^A^(^3[|i + ^S^fy-^ + C7) + BxiA^BWi + B^W2 + H^)^ 

-2p<p'^P-^UW - p{2p - ((2^J7V + U^){Bx + W) + A^^'^W) 

-P02P-2 (352^2^ ^ 35x^^ + W^) + iV„, 

Fh/ = Re{<l) + Aij + U) + {CiA^ + C2B'^)U + RBX + REdEX + N+. 



(4.16) 
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Substituting 
using the system 



and (|4.16|) into the system 
we obtain 



and computing the time derivative of {E, A, B) 



0, 
0, 
0, 



and 



0, 
0, 
0, 
0, 



and 



where 



L+iE)A + T + go,2iE) 
L_{E)T - 2G - 2A2(^)A + 



L+{E)Ur+k\E)Wi+g:ifl{E) 
L+{E)U2 + 2Wi + ?,K^[E)W2 + gi,2{E) 
L^iE)Wi - 3Ui - 2A\E)U2 - f2,i{E) 
L^{E)W2-U2- h,z{E) 



Ut = L^{E)W + Ru{E, A, B, U, W), 
-Wt = L+{E)U + Rw{E, A, B, U, W). 



Ru = Fu- {A^OeQ + B'^dE^ + A^OeUi + AB^dEU2)iC3AB + Re) 

- {2A@ + 3AUi + B'^U2){C^A'^B + C^B^ + Ra) 

- (25 A + 2ABU2){C(iA^ + CjAB^ + Rb), 

Rw = Fw + (ABdET + A'^BdEWi + B^dEW2)iC3AB + Re), 
+ (^r + 2ABWi){C^A'^B + C^B^ + Ra) 
+ (^r + A^Wi + 3B^W2){C6A^ + CyAB^ + Rb). 



(4.17) 



(4.18) 



(4.19) 



We shall now introduce two constrained spaces by 

Ll = {UeL\R): {4>,U)l2 = {x,U)l2=0}, (4.20) 
= {W£L\R): {dE(j),W)L2 = {ip,W)L2=0}. (4.21) 

Note that the orthogonal projections depend on E but we omit this dependence for the notational 
convenience. We can also define if|.(]R) as constrained spaces for any s > 0. In notations of Section 
3.1, we have used = H"^. 

The following two lemmas describe solutions of the systems (|4.17|) and (j4.18p . 

Lemma 4. There exists sufficiently small e > such that for all \E — E^\ < e, there exist a unique 
solution 0, A G //^(M) and T G //^(M) of the system (|4.17p . Moreover, these solutions are even in x, 
in E, and satisfy 

3C7>0: \\d%Q\\H2 + \\d%/^\\H2 + \\d%T\\H2 <C, q = 0,1,2. (4.22) 
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Proof. To solve the system (|4.17|) near E = E^, we recall that operators L^{E^) and L_{E^) are 
not invertible so that we shall set g2,o,9o,2 G and /i^i S thanks to the symplectic 

orthogonality. These constraints set up uniquely the coefficients Ci, C2, and C3, 

Ci=p2p+1 ,^ , , C2=p , , C3 = -2p — . (4.23 

{dE(l>,(p)L2 {dECpAjL^ {dECpAjL^ 

For E close to E'*, the existence and uniqueness of the solution 0, A € if^(R) and F G of the 

system (j4.17p follow from the gap between zero (or small) eigenvalues of L±{E) and the rest of the 
spectrum of L±{E) using the fixed point arguments. 

Indeed, for any fixed E such that \E — Ej,\ is sufficiently small, for any F G L?_(M), there exists a 
unique solution Gr, Ap G i7^(M) of the system 

L+{E)Qt + ^\E)T = -g2fl{E). L+{E)Ar + T = -go,2{E). (4.24) 

Then, for any two Fi,F2 G L^(R), we have 

L+{E) (en - GrJ = A^E) (F2 - Fi) , L+{E) (A^ - ArJ = F2 - Fi. 

So, by using the elliptic regularity, there is a positive constant C{E) such that 

IIGri - Gr2||H2 < C(E)|A(E)|2||F2 - FiH^a, ||Ar, - Ar^l^a < C{E)\\T2 - FiH^a. 

We recall that A^(S) — t- as ii^ — ?• ii^^, . From these estimates and by applying the fixed point arguments, 
we get the existence of F G (M) which solves the equation 

L_{E)T - 2Qr - 2k\E)l\r = -fi,i{E). (4.25) 

To estimate the H^-norm of the three solutions, we write equations (14.24p as 
L{E,)er = -g2,o{E) - {L+{E) - L+{E,)) Gr - A\E)r 

and 

L(S,)Ar = -go,2{E) - {L+{E) - L+{E,)) Ar - F. 

Because ||L+(£;)-L+(£'*)||loo =0{E-E^) and h?{E) = 0{E-E^) are sufficiently smah for | < 

e, the standard regularity theorem for elliptic equations implies that there exists an £'-independent 
constant C > such that 

||Gr|b2 <C(l + e||F||^2), ||Ar||i^2 < C (1 + ||F||j:^2) . (4.26) 

From these estimates and equation (j4.25p . we get for some constants C,C > 0, 

\\T\\h2 < C7(1 + ||Gr||//2 +e||Ar||H2) < C (1 + e||F||^2) . 

So, if \E — E^:\ < e is small enough, such that Ce < 1, there is an £^-independent constant C > such 
that ||F||j^2 < C. From this estimate and estimate ()4.26p . we also obtain the i/^-estimates of G and A. 

Because L±{E), g2fl{E), gQ^2{E), and fi^i{E) are all in E, we see that G, A, and F are all in 
E. By differentiating the system ()4.17p and using the same method as the one we just used, we also 
obtain the //^-estimates of 9^G, d'^A and 5gF with a = 1, 2. □ 
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Lemma 5. There exists sufficiently small e > such that for all \E — < e, there exist a unique 
solution Ui,U2 G and Wi,W2 G -ff^(M) of the system ( |^. jgp . Moreover, these solutions are odd 

in X, in E, and satisfy 



3C > 



3%Uj\\H2 + \\d%Wj\\H2<C, j = l,2, a = 0,1, 2. 



(4.27) 



Proof. To solve the system ()4.18p near E = E^, we set §3,0,51,2 G i^^(K) and /2,i, /o,3 G -^+(1^) thanks 
to the symplectic orthogonahty. These constraints set up uniquely the coefficients C4, C5, Cg, and C7, 

-C3(9i=;V,x)L2 +cri||x||i2 - 2p(x^0'^"'e)i2 -p(2p- l)(x^02^-V)L2 



C5 
C7 



C2\\X\ 



L2 



2 ^2p-l 



A) 



L2 



CillVlli^ +2p(2p + l)(^2^,/.2f-ie),.2 + :.p(2p + l)(2p- l)(^^</>2p-V>L2, 



1 

3 

+ p(2p-l)(^2^02p-2^2^^^^ 

The rest of the proof is similar to that of Lemma HI 



(4.28) 



□ 



To end this section, we estimate the error terms Rq, Re, Ra, and Rb of the system (j4.8 
Lemma 6. There is C > such that 

\Re,E,A,B\ < C ({A^ + B^f + {A + B){\\U\\l2 + \\W\\l2) + WWl^ + \\W\\l^ 



Re 
Re 



{(t),N_{A^ + U,Bx + W))L2 
{dE<p,N+{A,l: + U,Bx + W))l2 



Proof. We recall from (|iT3]) and (ji^ that 

{{dE(pA)ll+0{\\U\\L2 + \\W\\L2)) 

By using the symmetry properties of (p, ip, and as well as Lemmas H] and O we get 

{cl),N.{Ai; + U,Bx + W))L2 = {SecPA) M{E)AB + {4>,N.) l2, 

- (dEcj), N+{AiP + U,Bx + W)) L2 = {dEct>, 0) [CiA^ + Cs^^] - {OeP, 



(4.29) 



(4.30) 



(4.31) 



where Ci{El, C2iE)^and C-i{E) are defined in It then follows from (gSI), (ITOD . and 

(lOTT) that Re and satisfy (ltl9]) . 

The computation of the terms and -R^ can be done exactly the same way using ()4.5p . (|4.6p . and 
the above estimates on i?^; and Rq. □ 



4.2 Analysis of the remainder terms 



We shall now consider the remainder terms U and W satisfying the system ()4.19p . Let us denote 
Z = (i/, TV), R = {Ru, -Rw), and 



H{E) 



L_{E) 
-L+{E) 



21 



The system (|4.19p can be rewritten in the matrix-vector notations as 

OtZ = H{E)Z + n{E, A, B, Z). 



(4.32) 



Let Pc{E) be the projection operator associated to the complement of the four-dimensional subspace 
spanned by 

















X 







dE4> 












We shall single out quartic terms in {A,B) from the nonlinear term of the system (|4.32p . To do so, 
we expand 

R= ^ {ij{E)A'B^ + F{E,A,B,Z), (4.33) 

i+j=4 

where the functions %(S) = PciE)fij{E) G H^{R) can be computed explicitly near E = E^, 
whereas the function F{E, A, B, Z) satisfies the bounds 

3Cs>0: ||F||i/= (|^|^ + |S|5 + (|A| + |S|)||Z||iy. + ||Z||1^.), (4.34) 

for any s > (Recall here that H''^{M.) is a Banach algebra for any s > ^.) Without loss of generality, 
we can work for s = 1. 

Now, we mirror the decomposition ()4.33p and expand Z as 

Z= Zij{E)A'B^ + Z{E,A,B), (4.35) 

where the functions Zij{E) = Pc{E)zij{E) £ i/^(]R) can be computed explicitly near E = E^. The 
new variable satisfies 

dtZ = H{E)Z + -R{E,A,B,Z), (4.36) 

where the residual term R is computed from F similar to how R is computed from F. Therefore, the 
residual term satisfies the bound 

30 0: \\R\\Hi<c(^\Af + \B\^ + {\A\ + \B\)\\Z\\Hi + mliy (4.37) 

Because E depends on t, the spectral projections associated to the linearized operator H{E) are time 
dependent. Since E is close to i?*, we can fix the value E^, before writing the time evolution problem 
(j4.36p in the Duhamel form. In other words, we first rewrite (j4.36p as 

dtZ = P:{E,)H{E,)P,iE,)Z + Hr{E)Z + R(^, A, B, Z), (4.38) 

where Hr{E) = P*{E)H{E)Pc{E) - P*{E^)H{E^)Pc{E^) is a 2 x 2 matrix-valued function. Thanks 
to smoothness of this function in E and x, it enjoys the bound 

3C>0: \\Hr{E)\\ci<C\E-E^\. (4.39) 
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Then we use the Duhamel principle and write 
Z{t) = P,*(E,)e*^(^*)p,(ii;,)Z(0) 

+ P:{E,) /^(*-^)^(^*)Pe(£;,)f//«(£;(T))Z(r) + R(i?(T),yl(r),i?(r),Z(r))ldT. (4.40) 
Jo L ^ 

It follows from stability of the symmetric state <j){E) in Theorem [2] that the operator H{Ejf) has zero 
eigenvalue of algebraic multiplicity four and the rest of the spectrum is purely imaginary and bounded 
away from zero. Therefore, operator P*{E^)e^^^^*^ Pc{E^) forms a semi-group from ff*(M) to H^(K) 
for any s > and there is Cg > such that 

||P*(i?,)e*^(^*)p,(i?,)||H^^H^ < Cs. (4.41) 

Local existence and uniquness of solutions Z(t) of the integral equation ()4.40p follows for any t € 
[0,to]) where to > is sufficiently small for fixed-point arguments. The solution can be continued over 
the maximal existence interval using standard continuation methods. We shall now use Gronwall's 
inequality to control ||Z(t)||j:/i over t G [0, T], where T is bounded by (I4.1ip and {A,B) belong to the 
domain (|4.1U|) . 

Using (OT]) . (09]1 . (OOjl . and (IOT]l . we obtain for any t £ [0,T] that 

3O0: WmWHi < C||Z(0)||hi+C /Vw-^*ll|Z(r)|bidr + C /Vwil|Z(T)||Hidr 

JO Jo 

+ CT\ANf/^ + C [ 
Jo 

By Gronwall's inequality, we have 

sup ||z(t)||^i < C7 ( ||z(o)||^^i +r|Aiv|5/2^r sup ||z(t)||^i ) e^-/'o^(i-^(")-^*i+i^(")i)'^". 

i6[0,T] y te[0,T] J 

If A, T, and E are estimated by ()4.10p . ()4.1ip . and ()4.12p respectively, the last exponential term is 
bounded as |AiV| — t- 0. Elementary continuation arguments give that if ||Z(0)||j:/i < Co(AA^)^, then 
there is C > such that 

sup ||Z(t)||^i <C{ANf, te[0,T], 

t£[0,T] 

or, by virtue of equations (|4.10p and (j4.35p . there is C > such that 

sup ||Z(t)||j^i < C7(AiV)2, te[0,T], (4.42) 

tG[0,T] 

Bound ()4.42p provides the proof of the estimate ()4.13p . The estimate (j4.14p follows from (j4.29p and 
()4.42p . All together, the proof of Theorem [H] is complete. 



''*"Z(r)|||idr. 
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4.3 Conserved quantities 

To complete the proof of Theorem [3l we need to show that the trajectories of the system (j4.8p remain 
in the domain ([iTTU]) for t £ [0,T] and satisfy the estimates (jnT]) and (0712]) • 

Estimates (|4.12|) follow from the first two equations of the system (j4.8p in the domain (|4.1U|) under 
the estimate (j4.14p on the error terms and the estimate (|4.1ip on the maximal time T. Therefore, we 
shall only prove the estimates (|4.10p and (j4.1ip . To do so, we work with the last two equations of the 
system (j4.8p and employ the conserved quantities (j2.23p and ()2.24p . 

Expanded at the quadratic terms in (A, B), the conserved quantity for A/q becomes 
Mo = NsiE) + A'^\\ij\\l2 +B'^\\x\\l2+o(^iA^ + bY + {A + B){\\U\\l2 + \\W\\l2) + \\U\\l2 + \\W\\l2 

We note that the terms involving U and W are controlled by the bound (|4.13p to be of the higher 
order than the terms involving A and B in the domain (|4.10p . To simplify our notations, we shall then 
rewrite A/q simply as 

Afo = Ns{E) + A^Ml^ + B^xWh + ^(^' + B^f- (4-43) 

Remark 5. Computing the derivative of ^.43\ ) in time and using system |^.<g[ j up to the quadratic 
order, we obtain 

K{E)C,{E) + 2mE)\\l2 + K\E)\\x{E)\\l2) = 0, 
which is identically satisfied from system i2.12\) thanks to the identity 

2 I a2|UJ|2 



£2+A^||x||i2 = {ij,L_x)L^-{x,L+ij)L^ = {iJ,iL^-L+)x)L^ 

= 2p{i^,4>'"'x)L^ = -\K{E)C^{E). (4.44) 

Expanded at the quadratic terms in {A,B), the conserved quantity for Hq becomes 

Ho = Hs{E) + A^ [ {-^l + - (2p + 1)'^'^V-'') dx + B^ [ {xl + Fx' - ^^"x^) dx 

JR JR 

+0{A^ + B'^f, (4.45) 
Using the system (|2.12p and the normalization {iP,x)l'^ = 1) we obtain 

r (V;2. + yV''-(2p + l)W)dx = -A\E)-E\ml2, 

JR 

I {xl + Vx'-<p'Px^)dx = l-EWxWh- 

JR 

Using (j4.43p . we can further simplify ()4.45p to the form 

■Ho = Hs{E) + E{Ns{E) - Mo) - A^{E)A'^ + B"^ + 0{A^ + B^f. (4.46) 

We can now extend the conserved quantity Ho up to the quartic terms and write it abstractly as 
Ho = Hs{E) + E{Ns{E)-Mo)-A\E)A' + B^ 

+ ^Di{E)A'^ + D2{E)A'^B^ + ^D3{E)B'^ + 0{A^ + B^f, (4.47) 
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where Di, D2, and -D3 are some coefficients, which can be computed exphcitly. To avoid lengthy 
computations, we shah compute these coefficients from the derivative of (j4.47p in time and using the 
identity (|2.25|) and the system (|4.8|) up to the quartic order. This procedure yields two relations between 
three coefficients Di, D2, and D^, 

Di + A^D2 = ^CsiUWl^ + OeA^) + C^A^ - Ce, 
D2 + A^D; = ic3||x|li2+C5A2-C7. 

Since A^(i?*) = 0, coefficients Di{E^) and D2{E^,) are determined uniquely from this system. In 
particular, using (j2.10p . (j4.23p . and (|4.44p . we compute 

n^^^^ o 2/^ , ,-, {dE(l)*,(pf~^'^*)L2{i^*,<pfx*)L2 ^ 
Di{E^) = -2p{2p + l) ,^ , , , Q 

= -p{2p + i y ^ - q. (4.48) 



4.4 Analysis of dynamics as — )■ i?* 



Hamiltonian system ()4.8p equipped with conserved quantities (j4.43p and (|4.47p is integrable in the 
sense of the Liouville up to the error terms controlled by Theorem [H Using the conserved quantity 
(j4.43p and the assumption that N'^{E^) ^ 0, we can exclude the variable E near E = E^. Then, 
using the conserved quantity (14.47p . we can plot the trajectories of the system (I4.8P on the phase plane 
{A, B) and show the topological equivalence of the phase portraits of the system (|4.8p to those of the 
second-order system (ll.9p . 

By the assumptions of Theorem [3l we have X'{E^) < 0, N!,{E^) > 0, and Q < 0. As stated in 
Theorem \5\ these assumptions are satisfied for the symmetric potential (jl.2p in the limit of large 
separation of wells, that is, as s — t- oo. Hence the asymmetric states (p±{x;£) in the expansion (I3.14p 
exists for £ > E^. 

We now proceed with the phase plane analysis for the system ()4.8p . We denote AA^ = A/q — A^s(S*) 
and /S.H = T-Lq — Hs{E^:). Assuming that AA^ is small, we shall work in the domain ()4.10p and use the 
expansions ()4.43p and (I4.47P rewritten again as 

AN = N'^{E,){E - E,) + A^W^.Wl^ + B^\\xA\l2 + 0{{E - E,f, [E - E,){A^ + B^), {A^ + B^f) 
and 

AH = Hs{E)-Hs{E,) + E{Ns{E)-Ns{E,)- AN) 

- A^{E)A^ + B^ + ^Di{E)A'^ + D2{E)A^B^ + ^Ds{E)B^ + 0{{A^ + B^f), 

= -EAN + ^N'^{E,){E - E,f + \' {E,)\\iP4l-2{E - E,)A^ + B^ + ]^Di{E,)A^ + D2{E,)A^B'^ 

+ ]^D^{E,)B^ + 0{{E - E,f, {E - E,)\A^ + B^), {E - E,){A^ + B^f, {A^ + B^), 
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where we have used (|2.14|) and (|2.25|) . The first conserved quantity is useful to ehminate E in the 
domain ()4.10p by 

K{E,){E - E,) = AN- A^\\i;4l, " B^\\x*\\h + OiAN)\ (4.49) 
The second conserved quantity can now be written in the form 

G = i^I^^X'{EM*\\hA^ + B' + \f^A^ + F^A^B^ + ^FsB' + 0{ANf, (4.50) 

where 

(AiV)2 



G = AH + E^AN + 



2N',{E,y 



/ l|2 



F. . ^i^%|#,Z^ + Z..(^.), 



Using ([L5]), ([LT]), (l2l4l) . and ([QH]) . we obtain 

In the domain (14.101) . where = 0{AN) and B"^ = 0{ANy, G can be rewritten in the form, 

C = ^^X-iEMA^ . - ^^A' . OiANf. (4.51) 

This quantity is time-preserving for any t G [0, T]. Truncation of (j4.5ip at the error term OdAiVp) 
produces the second-order system ()1.9p . 



We can now fook at four different cases in the dynamics of the system produced by the conserved 
quantity (j4.5ip for any t E [0, T]. Recall that Ng{E^) > and Q < by our assumption. 



4.4.1 Case AiV > and cS > 

It follows from ()4.5ip that the critical point {A, B) = (0, 0) is a saddle point of G if AN > (recall 
that X'{E^:) < 0). The level set G = gives a curve on the phase plane {A,B) given by 



^ = + J^^A' + 0(ANf. (452) 

This curve contains the point (0,0) up to the terms of the order of (AA^)^. 
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If Q < and S > 0, the curve G = consists of two symmetric loops on the plane {A, B) that 
enclose the points (±74*,0), where G is minimal. An elementary computation shows that 



A' 



||2 




(4.53) 



QS 



It is easy to check that expansion (j4.53p follows from the stationary normal-form equation if 
[E — E^) is eliminated by the conserved quantity (j4.49p . 

The critical points (±^^,,0) are minima of G (center points). Therefore, they are surrounded by 
continuous families of periodic orbits on the phase plane {A^ B). Periodic orbits fill the domain enclosed 
by the two loops of the level G = 0. There are also periodic orbits for G > that surround all three 
critical points (ibA*,0) and (0,0). 

Thus, the truncated system ()4.8p approximates the dynamics near the unstable symmetric state 
<p{x;E) and the stable asymmetric states ip±{x]£). The phase portrait is topologically equivalent to 
the one on Figure [2] (top left). 

4.4.2 Case AA > and 5 < 

The critical point {A, B) = (0, 0) is again a saddle point of G if AA > 0. No other critical points of G 
for small A and B = Q exist if 5 < 0. Therefore, trajectories on the phase plane {A, B) near (0, 0) are 
all hyperbolic and they leave the neighborhood of (0, 0) in a finite time, except for the stable manifolds. 
The truncated system (j4.8p in this case approximates the dynamics near the unstable symmetric state 
<j){x;E). The phase portrait is topologically equivalent to the one on Figure [2] (top right). 

4.4.3 Case AA < and S >0 

It follows from ()4.5ip that the critical point (A, B) = (0, 0) is a minimum of G if AA < 0. No other 
critical points of G for small A and 5 = exist if 5 > 0. Therefore, (0,0) is a center point, which 
is surrounded by a continuous family of periodic orbits on the phase plane {A,B). The truncated 
system (j4.8p in this case approximates the dynamics near the stable symmetric state (f){x]E). The 
phase portrait is topologically equivalent to the one on Figure [2] (bottom left). 

4.4.4 Case AA < and 5 < 

The critical point (A, B) = (0, 0) is again a minimum of G if AA < 0. If 5 < 0, there exist two 
symmetric maxima of G at the points (±^*,0), where A^ is given by (j4.53p . The points (ibA*,0) 
are saddle points and the level curve ol G at A = A^^ and B = Q prescribes a pair of heteroclinic 
orbits connecting (=bA*, 0). The pair of heteroclinic orbits encloses a continuous family of closed curves 
surrounding the center point (0, 0) on the phase plane and corresponding to periodic orbits. The 
truncated system (j4.8p approximates the dynamics near the stable symmetric state 4'{x; E) and the 
unstable asymmetric states ip±{x]£). The phase portrait in this case is topologically equivalent to the 
one on Figure [2] (bottom right). 
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4.5 The end of the proof of Theorem [3] 



It follows from Section 4.4 that all nontrivial solutions of the system ()4.8p in the domain (j4.10p are 
topologically equivalent to the ones given by the second-order system (|1.9|) . 

The estimate (|4.11|) on the maximal time T > 0, during which the solutions remain in the domain 
(j4.10p . follows directly from the integration of the system (|1.9p over the time t E [0,r]. The other 
estimates of Theorem [3] follows from Theorem O The proof of Theorem [3] appears to be complete. 
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